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Analysis of gene expression from the Wolbachia genome 
of a filarial nematode supports both metabolic 
and defensive roles within the symbiosis 
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The a-proteo bacterium Wolbachia is probably the most prevalent, vertically transmitted symbiont on Earth. In contrast 
with its wide distribution in arthropods, Wolbachia is restricted to one family of animal-parasitic nematodes, the Oncho- 
cercidae. This includes filarial pathogens such as Onchocerca volvulus, the cause of human onchocerciasis, or river blindness. 
The symbiosis between filariae and Wolbachia is obligate, although the basis of this dependency is not fully understood. 
Previous studies suggested that Wolbachia may provision metabolites [e.g., haem, riboflavin, and nucleotides) and /or 
contribute to immune defense. Importantly, Wolbachia is restricted to somatic tissues in adult male worms, whereas females 
also harbor bacteria in the germline. We sought to characterize the nature of the symbiosis between Wolbachia and 0. ochengi, 
a bovine parasite representing the closest relative of 0. volvulus. First, we sequenced the complete genome of Wolbachia strain 
wOo, which revealed an inability to synthesize riboflavin de novo. Using RNA-seq, we also generated endobacterial 
transcriptomes from male soma and female germline. In the soma, transcripts for membrane transport and respiration 
were up-regulated, while the gonad exhibited enrichment for DNA replication and translation. The most abundant 
Wolbachia proteins, as determined by geLC-MS, included Iigands for mammalian Toll-like receptors. Enzymes involved in 
nucleotide synthesis were dominant among metabolism-related proteins, whereas the haem biosynthetic pathway was 
poorly represented. We conclude that Wolbachia may have a mitochondrion-like function in the soma, generating ATP for 
its host. Moreover, the abundance of immunogenic proteins in wOo suggests a role in diverting the immune system toward 
an ineffective antibacterial response. 



[Supplemental material is available for this article.] 

The a-proteobacterium Wolbachia is estimated to infect two-thirds 
of all arthropod species (Hilgenboecker et al. 2008) and thus is 
probably the most abundant vertically transmitted organism in 
the biosphere. However, it also has a much more restricted preva- 
lence in two families of nematode worms: about 20 species of an- 
imal-parasitic Onchocercidae (superfamily Filarioidea; the filariae) 
(Ferri et al. 2011) and at least two species of plant-parasitic Praty- 
lenchidae (Haegeman et al. 2009). The most salient characteristic 
of Wolbachia is its astonishing ability to induce a wide diversity of 
phenotypes in arthropod hosts, from reproductive manipulations, 
such as cytoplasmic incompatibility and male-killing, to mutual- 
istic roles, including pathogen protection and enhanced fecundity. 
Until recently, obligate dependencies on Wolbachia had not been 
recognized in arthropods, but a small number of species in dispa- 
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rate taxa require these symbionts for normal reproduction or de- 
velopment (Dedeine et al. 2001; Timmermans and Ellers 2009; 
Hosokawa et al. 2010; Miller et al. 2010). In contrast, obligate 
mutualism is widespread in the Onchocercidae, where antibiotic 
treatment induces retarded larval growth (Hoerauf et al. 1999), 
embryostasis in females (Bandi et al. 1999), and even death of adult 
filariae (Langworthy et al. 2000). Thus, this symbiotic relationship 
has been exploited as a drug target in the search for practicable 
adulticidal chemotherapy against human filarial pathogens, such 
as Onchocerca volvulus (the etiological agent of onchocerciasis, or 
river blindness) and Brugia malayi (one species responsible for 
lymphatic filariasis [LF]) (Hoerauf 2008). In adults of B. malayi, O. 
volvulus, and Onchocerca ochengi (a parasite of cattle representing 
the closest relative of O. volvulus) (Trees 1992; Morales-Hojas et al. 
2006), Wolbachia are found both in the somatic hypodermal cords 
that run along the length of the worms and in the germinal zones 
of the female gonad, where they are transmitted into the oocytes 
and maintained throughout larval development. However, they 
are completely absent from the reproductive system of adult males 
(Landmann et al. 2010). 
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Wolbachia are divided into about 10 supergroups on the basis 
of phylogenetic distance (Lo et al. 2007). Most symbionts in filarial 
nematodes are found in supergroups C, D (which both lack 
members in arthropods), or F (which also contains strains from 
bedbugs and termites, among other arthropods). To date, four 
Wolbachia genomes have been completed: two from group A 
(wMel, 1.3 Mb [Wu et al. 2004], and wRi, 1.5 Mb [Klasson et al. 
2009] from Drosophila melanogaster and Drosophila simulans, re- 
spectively), one from group B (wPip, 1.5 Mb [Klasson et al. 2008] 
from the mosquito Culex pipiens), and one from group D (wBm, 
1.1 Mb [Foster et al. 2005] from B. malayi). Comparative analysis of 
the wBm genome and that of B. malayi suggested that Wolbachia 
provisions the worm with haem, riboflavin, and flavin adenine di- 
nucleotide (FAD), since the B. malayi nuclear genome does not en- 
code complete pathways for these cofactors (Foster et al. 2005; 
Ghedin et al. 2007). In addition, and unusually for a member of the 
Rickettsiales, wBm retains the de novo nucleotide synthesis pathway 
and thus could also supply this resource during periods of high 
metabolic demand (Foster et al. 2005). Conversely, the almost 
complete lack of amino acid synthesis pathways in wBm clearly 
defined the nature of the mutualism in the reverse direction. More 
recently, an alternative or additional role for strain wOo (supergroup 
C) in O. ochengi has been proposed, as this symbiont stimulates an 
ineffective neutrophilic immune response in the bovine host (as 
also observed in human onchocerciasis) that apparently prevents 
the degranulation of eosinophils against adult worms (Hansen et al. 
2011). This may represent a form of defensive mutualism. 

Here, we present the most comprehensive molecular analysis of 
a single Wolbachia strain performed to date. First, we sequenced the 
complete genome of wOo and determined that it is the smallest yet 
described for the taxon. Second, we obtained complete transcriptomes 
from Wolbachia for the first time and compared the global gene ex- 
pression between somatic and gonad tissue. Finally, we identified the 
most abundant protein species in wOo from adult female worms. 
Our data reveal important dissimilarities between wBm and wOo in 
terms of metabolic capability, the density of insertion sequences 
(ISs), and the repertoire of repeat-motif containing proteins. More- 
over, our gene expression analyses indicate potential roles for wOo in 
both energy production in the somatic tissue and modulation of the 
mammalian immune response but fail to provide strong support for 
the provisioning of vitamins or cofactors by this strain. 

Results 

The genome of Wolbachia strain wOo is highly degraded 

At 0.96 Mb, the wOo genome is the most reduced Wolbachia ge- 
nome sequenced to date (11% smaller than that of wBm), and it 
also exhibits the lowest G+C content and gene density for any 
Wolbachia genome (Fig. 1; Supplemental Tables SI, S2). Although 
wOo shares several regions of microsynteny with wBm, gross 
chromosomal organization is markedly different between these 
strains, while conservation of gene order is reduced to the scale of 
only two to five genes between wOo and the Wolbachia genomes of 
insects (Supplemental Fig. SI). In addition, wOo exhibits the most 
extreme GC skew of any Wolbachia genome (Supplemental Fig. 
S2), which may reflect ancient gene order stability and, conse- 
quently, the accumulation of mutations on each DNA strand 
(Rocha 2004). The features that have been lost during reductive 
evolution of the wOo genome include transposable elements. 
Indeed, we only detected six IS copies (Supplemental Table S3), 
accounting for <0.5% of the genome, and no group II introns. This 
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Figure 1. The chromosomal structure of wOo. Circles are numbered 
sequentially from the perimeter to the center. Circles 1 and 2 represent 
protein-coding genes on the positive and negative strands, respectively, 
with intact CDS colored blue and pseudogenes in red; the third circle 
shows structural RNA genes (rRNA, purple; sRNA, green; and tRNA, pink); 
circles 4-7 are heat-maps of transcript counts in male soma (4, 5) and 
female gonad (6, 7), with high expression colored red and low expression 
in blue; circle 8 is a plot of differential expression (DE) between male soma 
(outer profile) and female gonad (inner profile), where red coloration in- 
dicates statistically significant DE at P < 0.05; and circle 9 represents the 
percentage of GC skew [(G - Q/(G + C)] for each DNA strand. Note that 
expression from tRNA and rRNA genes was forced to baseline to facilitate 
visualization of mRNA expression. 



is in marked contrast with all other Wolbachia genomes, in which 
ISs (in wBm and insect strains) and group II introns (in insect 
strains only) are found in high abundance (Supplemental Table S4; 
Cerveau et al. 2011; Leclercq et al. 2011). 

To define the unique set of protein-coding genes that are 
absent in strain wOo but present in all of the other completed 
Wolbachia genomes, we performed an orthologous cluster analysis. 
Genes that were apparently retained in wOo were then further 
characterized as intact or pseudogenized. In total, 88 genes have 
been uniquely lost or pseudogenized in the wOo genome, of which 
29 encode hypothetical proteins (Fig. 2). The remaining genes with 
a robust or putative annotation are predicted to be involved in 
DNA replication and repair, membrane transport, potential in- 
teractions with eukaryotes (repeat-motif proteins), and enzymes 
from cof actor synthesis and secondary metabolism pathways. Al- 
most every component of the homologous recombination ma- 
chinery is absent or pseudogenized in the wOo genome, including 
all members of the RecFOR repair pathway and the RuvABC 
complex that resolves Holliday junctions, and proteins conserved 
in almost all bacteria that could potentially complement these 
losses (RecA and RecG). The repertoire of membrane transporters is 
also somewhat depleted in the wOo genome. However, all subunits 
of the F D Fi ATPase and the ABC transporters that mediate the uptake 
of zinc (ZnuABC), ferric iron (AfuABC), and phosphate (PstSCAB) 
have been retained in wOo, alongside a Na+/H+ anti-porter that is 
absent in wBm. Moreover, we determined that the suite of genes 
associated with peptidoglycan synthesis is identical between wOo 
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Figure 2. Gene loss from the wOo genome. Orthologous cluster analysis identified 88 genes that are present in wBm, wMel, wPip, and wRi but that are 
lost or pseudogenized in wOo. These genes were manually classified into functional categories using the National Center for Biotechnology Information 
Conserved Domain Database. 



and wBm (Henrichfreise et al. 2009), with the important exception 
of the bifunctional N-succinyldiaminopimelate-aminotransferase/ 
acetylornithine transaminase gene, argD, which is pseudogenized in 
wOo (Supplemental Table S5). 

As cof actor metabolism is of particular interest in the context 
of the nutritional provisioning hypothesis, we examined the 
metabolic capabilities of wOo in relation to those of wBm. Notably, 
in contrast with wBm (Foster et al. 2005), most enzymes in the 
riboflavin pathway (including the final component, FAD synthase) 
have been lost or pseudogenized in wOo, with only the first gene in 
the pathway (ribA) remaining intact. Furthermore, the incomplete 
pathways for biotin, pantothenate, and coenzyme-A biosynthesis 
in wBm have been almost entirely lost from the wOo genome, and 
isoprenoid synthesis (secondary metabolism) (Fig. 2) is also se- 
verely restricted. However, all seven components of the haem 
pathway found in the other Wolbachia genomes are present in 
wOo. We note that in addition to the canonical enzymes of this 
pathway, all of the completed Wolbachia genomes contain a gene 
for HemY, a membrane-bound protein that can oxidize proto- 
porphyrinogen IX in Bacillus subtilis (Hansson and Hederstedt 
1994). Thus, hemY may complement protoporphyrinogen-IX ox- 
idase (hemG), which is absent in Wolbachia and most other 
a-proteobacteria (Foster et al. 2005). The wOo genome also contains 
the six cytochrome c maturation (can) genes present in wBm and 
the insect strains, which are responsible for the covalent attachment 
of haem to apocytochromes; whereas the gene for bacterioferritin 
(Bfr), a haem-storage protein, has been pseudogenized. We identi- 
fied, in contrast with the overall pattern of advanced genomic 
degradation in wOo, a downstream folate synthesis pathway in this 
genome, which is also present in wRi, but not in wBm, wMel, and 
wPip. As folate plays a key role in the methylation of DNA, these 
differences suggest that Wolbachia strains vary in their dependency 
on scavenging from the host for this essential cofactor. 

To determine whether the genomes of wOo and wBm lack 
common features that are conserved among all of the complete 
genomes from insect strains, we defined a repertoire of orthologous 



clusters that are shared among wMel, wPip, and wRi only (Supple- 
mental Fig. S3). This analysis confirmed the current paradigm that 
the major genomic distinctions between insect and filarial strains of 
Wolbachia can be attributed to prophage-related genes and those 
encoding repeat-motif proteins (especially ankyrin domains). We 
also calculated that the core set of orthologs shared between all of 
the completed Wolbachia genomes, including wOo, totals 5 75 genes 
(Supplemental Table S6). A comparison between this core set and 
those of Holman et al. (2009) (based on computational rankings of 
the essentiality of genes in the wBm genome by two different 
methods) displayed good concordance, as only three genes with 
high essentiality scores were absent from a Wolbachia genome 
(Supplemental Table S7). Notably, all of the genes encoding the 
type-IV secretion system are also intact in wOo (Supplemental Table 
S8) and conserved across all other Wolbachia genomes (Pichon et al. 
2009). 

The Wolbachia transcriptome reveals a constitutive 
heat-shock response 

To characterize the Wolbachia transcriptome, we performed RNA- 
seq on both whole adult male worms and gonads microdissected 
from adult female worms. Although only —5% of reads mapped to 
the wOo genome, almost 98% of predicted gene models (whether 
intact or pseudogenized) were classified as "expressed" (i.e., a 
minimum of 10 cumulative reads were detected across four in- 
dependent biological samples). However, gene models containing 
an intact CDS were expressed at a significantly higher level than 
were pseudogenes (P< 0.001) (Supplemental Fig. S4). To determine 
the distribution of gene function across the transcriptome, we 
classified protein-coding genes in the wOo genome using a com- 
bination of databases (Fig. 3). This analysis revealed that protein 
metabolism (including translation) was the most highly expressed 
functional category, followed by RNA metabolism (including 
transcription) and stress response. The latter subsystem was dom- 
inated by proteins that mediate redox reactions via thiol groups, 



Genome Research 2469 



www.genome.org 



Darby et al. 



□ Amino Acids and Derivatives 
□carbohydrates 

□Cell Division and Cell Cycle 
■cell Wall and Capsule 

□Cofactors, Vitamins, Prosthetic Groups, Pigments 
HDNA Metabolism 

□ Fatty Acids, Lipids, and Isoprenoids 
□Membrane Transport 
□Miscellaneous 

■ Nitrogen Metabolism 

□ Nucleosides and Nucleotides 



□Phages, Prophages, Transposable elements, Plasmids 

■Phosphorus Metabolism 

□Potassium Metabolism 

□ Protein Metabolism 

□Regulation and Cell signaling 

■Respiration 

□RNA Metabolism 

□stress Response 

■(Jncharacterised 

□virulence, Disease and Defense 




Figure 3. Classification of the wOo transcriptome. The wOo genome was categorized using a combination of subsystems from the SEED database, 
additional pathway descriptions from the Kyoto Encyclopaedia of Genes and Genomes, and manual curation based on the National Center for Bio- 
technology Information Conserved Domain Database. The number of transcript reads per functional category is shown in each slice of the pie-chart and 
represents the mean (i.e., reads are normalized to the total number of genes per subsystem). Exploded slices show (in clockwise order from 1 2 o'clock) 
"cofactors, vitamins, prosthetic groups, pigments"; "nucleosides and nucleotides"; and "protein metabolism." 



including thioredoxins, peroxiredoxins, glutaredoxins, and glu- 
tathione S-transf erases. Of the two functional categories that have 
been proposed to support nutritional provisioning in the symbi- 
osis with filariae, nucleoside and nucleotide metabolism was well 
represented, whereas pathways for the biosynthesis of cofactors 
and vitamins were expressed at a relatively low level. Remarkably, 
of the highly abundant transcripts ranked within the top 20, five 
(groEL, groES, clpB, dnaK, and Ion) are for chaperones or proteases 
associated with the heat-shock response (Supplemental Table S6), 
indicating constitutive activation of this pathway. 

Differential expression analysis supports limited regulation 
of growth rate and membrane transport between soma 
and germline 

To illuminate the potentially distinct roles for Wolbachia in the 
adult female germline and the somatic hypodermal cords of O. 
ochengi, we compared global transcription in wOo between these 



tissues. In accordance with the paucity of transcriptional regula- 
tors in the wOo genome, >96% of genes displayed a similar level of 
expression between the two anatomical sites (Fig. 1). However, we 
identified 26 genes that were differentially expressed at a signifi- 
cance threshold of P < 0.05 and a false-discovery rate (FDR) of 5%, 
encompassing fold-changes (FCs) within the range of 1.9-7.0. It 
should be noted that because our analysis was restricted to two 
biological replicates per tissue, it did not constitute a definitive 
experiment, but rather proof-of-principle that RNA-seq can be 
used in this context with obligate intracellular organisms despite 
a vast excess of host transcripts. In addition, the magnitude of 
these transcriptional differences was relatively modest, although it 
was comparable to that reported for Wolbachia ankyrin gene ex- 
pression between male and female gonads in Drosophila spp. as 
measured by quantitative real-time PCR (qRT-PCR) (Papafotiou 
et al. 2011). 

The pattern of differential expression was predominantly as- 
sociated with three functional categories: translation machinery, 
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DNA replication, and membrane transport (Table 1; Supplemental 
Table S9). In the germline, elevated expression of genes for elon- 
gation factor Tu (both paralogs), three ribosomal proteins, and the 
5S rRNA constituted a clear signal of increased protein synthesis. 
Of these, rpsU encoding ribosomal protein S21 exhibited the 
greatest FC of any transcript. Although rpsU is part of the macro- 
molecular synthesis operon in the Enterobacteriaceae, this is not the 
case in Wolbachia and many other eubacteria. For instance, ex- 
pression of an alternative rpsU gene cluster in cyanobacteria is 
regulated by temperature, such that the levels of protein S21 be- 
come substoichiometric at 38°C, despite optimal cell growth (Sato 
et al. 1997). This suggests that S21 is not essential for ribosomal 
function but may have a regulatory role in protein biosynthesis. 
Two nonstructural components with ancillary roles in translation 
were also up-regulated in the gonad: polypeptide deformylase 
(Mazel et al. 1994) and a predicted S-adenosyl-L-methionine-de- 
pendent methyltransferase (Pleshe et al. 2005; Savic et al. 2009). 
Additionally, evidence for a higher rate of DNA replication in this 
tissue relative to the soma was apparent in the increased expression 
of uracil-DNA glycosylase, a key enzyme of the base excision sys- 
tem (Holmquist 1998), and a nucleoid DNA-binding protein of the 
HU family, which exhibits peak expression during the exponential 
phase of growth in Escherichia coli (Luijsterburg et al. 2006). 
Moreover, the up-regulated gene annotated as a "cardiolipin syn- 
thase" has a conserved domain with closer homology with Nuc, an 
EDTA-resistant nuclease that digests nucleic acids without se- 
quence specificity (Pohlman et al. 1993). Since Nuc has a signal 
peptide, it could be involved in the recycling of nucleotides in the 
host germline to facilitate symbiont cell division. Taken together, 
the repertoire of genes displaying elevated expression in the gonad 
is fully consistent with the current understanding of Wolbachia 
biology in filarial nematodes, as wBm divides rapidly during the 
mitotic proliferation of host oogoniae, and this precedes fertiliza- 
tion of the mature oocytes (Landmann et al. 2010). 

Membrane transport, particularly in the periplasmic com- 
partment, was also differentially regulated between the gonad and 
soma. Whereas two genes involved in rate control of preprotein 
translocation via the general secretion (Sec) pathway were up- 
regulated in the germline, SecF (Duong and Wickner 1997) and 
signal peptidase I (Paetzel et al. 2002), the somatic wOo tran- 
scriptome was enriched for metal ion transporters and accessory 
proteins from the respiratory chain (Table 1). These included 



a cation diffusion facilitator (a family of secondary divalent cation 
filters associated with resistance to heavy metals), which can 
transport zinc, iron, cobalt, nickel, and possibly manganese, but 
not calcium, magnesium, or copper (Nies 2003). Oxidative phos- 
phorylation has a particularly high demand for heavy metal ions, 
and indeed, a cytochrome c oxidase subunit together with an en- 
zyme that may facilitate cytochrome maturation (thiol-disulfide 
isomerase) (Sambongi and Ferguson 1994) was up-regulated in the 
soma. Furthermore, two additional components of the respiratory 
chain exhibited an FC of more than 1.8-fold in this tissue, al- 
though these differences narrowly failed our FDR cut-off (Supple- 
mental Table S6). 

Remarkably, two sRNA molecules were among the most dif- 
ferentially expressed transcripts. The up-regulation of a transcrip- 
tional regulator in the germline, 6S RNA, was unexpected as it at- 
tains maximum concentrations in E. coli cells as they enter late 
stationary phase. Alternatively, as 6S RNA appears to play a critical 
role in translational adaptation under conditions of nutrient stress 
(for which a primary signal is nucleotide availability ) (Wassarman 
2007), increased production in the gonad may reflect the imper- 
ative for Wolbachia to synchronize its replication rate with that of 
the host. In accordance with this hypothesis, 6S RNA is required for 
optimal intracellular replication of Legionella pneumophila (Faucher 
et al. 2010). The elevated expression of RNase P in the soma was 
also surprising, considering its canonical function in the catalysis 
of precursor tRNA maturation (Kazantsev and Pace 2006). How- 
ever, RNase P is a remarkably pleiotropic ribozyme that can regu- 
late the expression of the Sec pathway (Li and Altman 2003) and 
cleave transient secondary structures in riboswitches (Altman et al. 
2005). Thus, the loss of three transcriptional regulator proteins in 
wOo (Fig. 2) may have resulted in a greater dependency on sRNA 
mediators in this role. 



The profile of abundant proteins includes Iigands 
for mammalian Toll-like receptors 

To determine the degree of concordance between transcript 
abundance and protein expression and to identify protein prod- 
ucts from hypothetical genes, we conducted a proteomic anal- 
ysis of Wolbachia-emiched material from whole adult female 
O. ochengi. Unique peptides from a total of 122 wOo proteins were 
robustly identified (P < 0.01, 1% FDR), representing 18.9% of 



Table 1. Up-regulation of Wolbachia periplasmic protein-encoding genes in the Onchocerca ochengi soma 



Locus tag 


Annotation 


Fold-change 


Putative function 3 


wOo08460 


Hypothetical protein 


3.06 


TrbC/VirB2 pilin subunit of type IV 
secretion system 


wOo08840 


Major facilitator super-family 
permease 


2.66 


Secondary active transporter of inorganic 
ions or small organic molecules 


wOoOl 770 


Hypothetical protein 


2.46 


Multiple resistance and pH regulation protein 
MrpF; Na + /H + or K + /H + anti-porter 








wOo07860 


Hypothetical protein 


2.22 


Predicted transmembrane protein with signal 
peptide; function unknown 


wOo08890 


Cytochrome c oxidase, 
subunit 3 


2.18 


Stabilizing component of complex IV of the 
respiratory chain 


wOo03390 


Thiol-disulfide isomerase 


2.16 


Thioredoxin family; attachment of haem 
to apocytochromes 


wOo03990 


Co Zn Cd efflux system 
component 


2.13 


Cation diffusion facilitator family; efflux of heavy 
metal ions (ionic radii, 76-97 pm) 
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predicted proteins encoded by an intact CDS (Supplemental Tables 
S10, Sll). There was a statistically significant relationship between 
transcript abundance and the detection of corresponding peptides, 
in that the expression level of genes displaying peptide evidence 
was significantly higher than for genes lacking detectable peptides 
(P< 0.001) (Supplemental Fig. S5). However, we found no evidence 
of protein production from pseudogenes, even where they were 
highly expressed at the transcript level. We used the number of 
unique peptides detected as an approximate measure of abun- 
dance (Supplemental Table Sll), and this revealed that GroEL 
completely dominated the wOo proteome, with almost twice as 
many unique hits than any other protein. Thus, in common with 
other obligate intracellular bacteria, Wolbachia GroEL displaces 
elongation factor Tu as the most highly expressed gene product at 
the protein level. Importantly, this phenomenon has been de- 
scribed both for other members of the Rickettsiales (Hajem et al. 
2009) and more distantly related Proteobacteria (Aksoy 1995; 
Baumann et al. 1996), irrespective of the nature of the symbiosis. 
In both Buchnera and the unrelated Blattabacterium, GroEL appears 
to be under positive selection, and it has been hypothesized that its 
substrates may have greatly expanded in these obligate in- 
tracellular bacteria (Fares et al. 2005). Contrarily, we found no 
evidence for positive selection in Wolbachia GroEL (Supplemental 
Table S12). 

The dominant metabolism-related subunits in the wOo pro- 
teome were from succinyl-CoA synthetase and ATP synthase, 
which generate purine nucleoside triphosphates from the tri- 
carboxylic acid cycle and the respiratory chain, respectively. In 
addition, succinyl-CoA is an essential precursor for the biogenesis 
of haem. Intense activity at the earlier stages of nucleotide syn- 
thesis was reflected by the abundance of glycine-serine hydroxy- 
methyltransf erase, the major driver of the one-carbon pool by 
folate, and of the gluconeogenic enzyme, fructose-bisphosphate 
aldolase, that produces nucleoside backbones via the nonoxidative 
pentose phosphate pathway. Overall, and in agreement with 
metabolic reconstruction in wBm (Foster et al. 2005), the wOo 
proteome suggests that amino acids are the principal carbon 
source, and evidence for the scavenging of sulfur from methionine 
and cysteine was also apparent (Supplemental Table SI 3). In con- 
trast with nucleotide and amino acid metabolism, biogenesis of 
cofactors and vitamins was poorly represented, with few unique 
peptides detected for enzymes in the folate (FolKP) and haem 
(HemA, HemC) synthesis pathways. 

To perform a comparison of the predominant protein com- 
position in wOo and wBm, we ranked Wolbachia proteins by the 
number of unique peptides detected in adult female worms of B. 
malayi (Bennuru et al. 2011) and O. ochengi (Supplemental Table 
Sll). Remarkably, this rank ordering showed substantial differ- 
ences, even on the simple measure of presence in the top 20; al- 
though GroEL, DnaK, a single ortholog of Wolbachia surface pro- 
tein (WSP; wOo07430/wBm0432), and a putative outer membrane 
protein (wOo08110/wBm0010) predicted to be a porin (Supple- 
mental Table SI 4) clearly dominated the proteome in both strains. 
A key conserved feature of wOo and wBm is the abundance of three 
proteins (GroEL, WSP, and peptidoglycan-associated lipoprotein 
[PAL]) that are known to stimulate potent immune responses in 
the mammalian host. Indeed, GroEL and its eukaryotic homolog 
(HSP60) from numerous organisms have been demonstrated to in- 
duce innate immune signaling and autoimmune reactions (Maguire 
et al. 2002). For instance, GroEL is the immunodominant antigen in 
Rickettsia conorii (Renesto et al. 2005), and Porphyromonas gingivalis 
GroEL is a confirmed ligand for mammalian Toll-like receptor (TLR)-2 



and -4 (Argueta et al. 2006). Although the potential for TLR signaling 
by Wolbachia GroEL has not been investigated in filarial diseases, LF 
patients exhibiting chronic pathology showed higher IgGl antibody 
levels to wBm GroEL than did exposed but asymptomatic groups 
(Suba et al. 2007). However, GroEL may also contribute to the T-cell 
anergy observed in chronic LF via elevated production of the anti- 
inflammatory cytokines interleukin (IL)-10 and transforming growth 
factor-^ (Shiny et al. 2011). A much larger body of data exists per- 
taining to the immunogenicity of WSP. This molecule is a ligand for 
TLR-2 and -4 (Brattig et al. 2004) and stimulates production of the 
pro-inflammatory cytokine IL-8 in neutrophils (Bazzocchi et al. 2003) 
while inhibiting apoptosis in these cells (Bazzocchi et al. 2007). 
Moreover, in common with GroEL, antibody levels against WSP are 
conelated with the disease manifestations of LF (Punkosdy et al. 
2003). Finally, PAL has been recently identified as a major inducer of 
innate pro-inflammatory immune responses in filarial infections. 
A synthetic diacylated analog of PAL from wBm stimulated TLR-2 and 
-6, leading to secretion of tumor necrosis factor-a and IL-8 from 
mammalian cells in vitro and recruitment of neutrophils into the 
cornea in the murine keratitis model of River Blindness (Turner et al. 
2009). 

Diacylation, but not triacylation, of lipoproteins was pre- 
dicted in wBm as its genome encodes genes for prolipoprotein 
diacylglyceryl transferase and lipoprotein signal peptidase, but not 
apolipoprotein aminoacyl transferase (Turner et al. 2009). We de- 
termined that wOo has maintained an equivalent repertoire of li- 
poprotein biosynthesis genes, and thus sought to identify the 
complete coding potential for lipoproteins using a variety of in 
silico prediction tools. Intriguingly, despite the reduced genome 
size of wOo relative to wBm, eight putative lipoproteins were 
identified in the former and only five in the latter, in part due to 
lack of sequence conservation between paralogs of the type-IV 
secretion system component, v/rB6 (Supplemental Table SI 5). 
Furthermore, one of the predicted lipoproteins common to both 
strains, a putative SCOl/SenC family protein, was moderately 
abundant in the wOo proteome and could constitute an additional 
trigger for TLR2/6 signaling. This group of proteins represents 
thiol-disulfide oxidoreductases and plays an important role in the 
assembly of the Cu A redox center of cytochrome c oxidase in a di- 
verse range of organisms (Banci et al. 2011). 

Discussion 

Three principal hypotheses have been proposed to account for the 
symbiotic relationship between Wolbachia and the Onchocercidae, 
which are not mutually exclusive. First, Wolbachia may provide 
metabolites to the worm host that either are unavailable from al- 
ternative sources or could become limiting at certain stages of the 
lifecycle, thus functioning as a nutritional mutualist (Foster et al. 
2005). Second, the endobacteria could indirectly interact with the 
immune system of the definitive (and perhaps the intermediate) 
host(s) of these filarial parasites, contributing to immune evasion 
as a defensive mutualist (Hansen et al. 2011). Finally, the obligate 
dependency of some onchocercids on Wolbachia may not be a 
mutualism at all, but rather an extreme form of reproductive par- 
asitism that has successfully hijacked normal cellular development 
in the embryo (Landmann et al. 2011). These hypotheses must 
contend with two recent and interrelated developments in the 
field: the demonstration of secondary loss of Wolbachia in several 
members of the Onchocercidae (McNulty et al. 2010), and the 
revelation that Wolbachia is currently present in <40% of this 
taxon (Ferri et al. 2011). 
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The main metabolic candidates that have been proposed to be 
provisioned by wBm are haem, nucleotides, riboflavin, and FAD. 
The latter two cofactors can be discounted in the case of wOo, since 
the necessary pathway has been pseudogenized. In contrast, a ca- 
pacity for de novo haem biosynthesis appears to be universally 
conserved in Wolbachia but is absent (or at least incomplete) both 
in free-living nematodes such as Caenorhabditis spp. (Rao et al. 
2005) and in parasitic species, including B. malayi (Ghedin et al. 
2007). Thus, Wolbachia may have become established in the 
Onchocercidae as a ready source of this essential compound. 
However, it is important to emphasize that haem autotrophy is 
almost universal in bacteria (with the exception of some special- 
ized species that reside in mucosal tissue) (Smalley et al. 2011) and 
appears to be conserved across the Rickettsiales, irrespective of the 
phenotypic relationship of the bacterium to its host. In addition, 
ingestion of red blood cells in the filaria Litomosoides sigmodontis 
has been recorded (Attout et al. 2005), while the intestine of adult 
O. volvulus (an organ lacking bacteria) contains hemosiderin 
(George et al. 1985), suggesting that oral uptake of erythrocytes 
also occurs in this species. Despite these considerations, several 
lines of evidence indicate a key role for Wolbachia in the mainte- 
nance of iron homeostasis (as opposed to simply provisioning). In 
D. melanogaster, wMel was shown to confer a fecundity benefit on 
flies reared on both ecologically relevant, iron-restricted diets and 
on artificial, iron-overloaded diets (Brownlie et al. 2009). More- 
over, in D. simulans reared under high-iron conditions, flies 
infected with wRi absorbed more iron than did uninfected flies, yet 
the up-regulation of host ferritin was lower than that observed in 
uninfected flies (Kremer et al. 2009). In this system, increased ex- 
pression of Bfr by wRi apparently compensated for the host re- 
sponse, although this mechanism cannot apply to wOo, which has 
lost the Bfr gene. However, up-regulation of a cation diffusion fa- 
cilitator by wOo in host soma could represent an alternative 
symbiont strategy to control iron metabolism, in which the bac- 
teriophorous vacuole could serve as a buffer zone mediating ex- 
change between host ferroproteins and iron cations. 

The importance of nucleotide metabolism was strongly sup- 
ported both by the retention of key pathways in the wOo genome 
and the high maintenance of expression of relevant genes. In 
contrast with haem biosynthesis, the capacity to synthesize nu- 
cleotides de novo has been lost by pathogenic rickettsiae (Foster 
et al. 2005), suggesting a divergence that is relevant to the evolu- 
tionary development of parasitism or mutualism. However, the 
high profile of ATP synthase in the transcriptome and proteome, 
together with the elevated expression of respiratory chain com- 
ponents by wOo in the hypodermal cords, suggests that the prin- 
cipal molecule delivered to the worm host is ATP for energy rather 
than other nucleoside triphosphates for nucleic acid synthesis. 
Indeed, in accordance with this hypothesis, one of the early signs 
of morbidity in antibiotic-treated adult worms is a reduction in 
motility (Gilbert et al. 2005). Moreover, ultrastructural observa- 
tions of the hypodermal cords in adult O. ochengi indicate that 
Wolbachia density in infected cells greatly exceeds that of mito- 
chondria (Gilbert et al. 2005), perhaps reflecting a functional shift 
from the organelle to the symbiont for energy generation. In fur- 
ther support of this proposal, tetracycline treatment induced up- 
regulation of mitochondrion-encoded respiratory chain subunits 
in L. sigmodontis (Strubing et al. 2010), which may represent an 
attempt by the host to compensate for reduced ATP availability. 

Whether Wolbachia has a mutualistic role in the female 
germline is less clear. It has been suggested that wBm supplies 
nucleotides to the worm host during the period of the lifecycle 



when the demand for these molecules is greatest; i.e., during em- 
bryogenesis (Foster et al. 2005). Our gene expression data for wOo 
within gonad tissue, which revealed the up-regulation of 6S RNA 
and a secreted nuclease, are more compatible with intense com- 
petition between worm and symbiont for nucleotides during oo- 
genesis. However, following fertilization of the oocyte, the pres- 
sure on this resource may decrease, since the total number of 
bacteria per individual host remains relatively constant during 
embryonic and larval development, up to the third molt (McGarry 
et al. 2004). This raises the question of how the infection in the 
oogoniae is maintained when it could have a deleterious effect on 
the host. In the wasp genus Asobara, various levels of dependency 
on Wolbachia exist, such that all known wild populations of Aso- 
bara tabida require strain wAtab3 to complete normal oogenesis 
(Kremer et al. 2010). This phenotype appears to be the result of 
a critical interaction between Wolbachia and an apoptosis check- 
point in the nurse cells of the ovarioles (Pannebakker et al. 2007) 
because removal of the symbiont triggers widespread host cell 
death in this tissue. As a minority of antibiotic-cured lines of A. 
tabida can produce very small numbers of viable offspring (Kremer 
et al. 2010), this obligate dependency is probably evolutionarily 
recent. Remarkably, the severity of the ovarian phenotype in 
A. tabida is correlated with the expression of host genes involved in 
iron homeostasis and oxidative stress (Kremer et al. 2010). Recent 
studies in B. malayi have demonstrated that apoptosis in the fe- 
male reproductive tract also follows symbiont depletion by anti- 
biotic treatment (Landmann et al. 2011). 

Taken together, these observations suggest that infection of 
the female germline by Wolbachia can foster dependency relatively 
rapidly without any obvious fitness benefit to the host, thus rep- 
resenting an "obligate parasitism." This could explain the exis- 
tence of filarial species that lack symbionts in the hypodermal 
cords (Ferri et al. 2011), which may be at an early stage of contact 
with Wolbachia, and of hosts that were once infected with Wol- 
bachia but have become aposymbiotic (such as Onchocerca flexuosa) 
(McNulty et al. 2010), which could reflect the evolution of coun- 
termeasures to regain control of the apoptotic cascade. In addition, 
if Wolbachia in filarial nematodes is fundamentally a parasitic in- 
fection that has secondarily acquired mutualistic traits, the hy- 
pothesis that lateral gene transfers of symbiont DNA into the host 
genome could constitute the capture of useful functions is unlikely 
to be correct (Fenn and Blaxter 2007). These considerations lead us 
to conclude that the mutualistic phenotype in filariae arose as 
a result of Wolbachia colonizing the somatic hypodermal cords. 
The most probable fitness benefit was an increase in the re- 
productive lifespan of the worms, since females of O. volvulus (and 
almost certainly O. ochengi) remain fecund for >10 yr, whereas 
O. flexuosa is thought to only live for ~1 yr (Plenge-Bonig et al. 
1995). As reactive oxygen species are one of the principal drivers of 
aging (Buttemer et al. 2010), we propose that the ability of Wol- 
bachia to perform aerobic respiration and to metabolize iron while 
maintaining an elevated oxidative stress response are key mecha- 
nisms that promote infection of somatic tissues in filariae. 

The longevity of O. volvulus and O. ochengi is also dependent 
on an ability to withstand the local immune response of the de- 
finitive host while enclosed within a fibrous nodule. In both of 
these species, Wolbachia is responsible for the recruitment of an 
ineffective neutrophilia around the adult worms (Brattig et al. 
2001; Nfon et al. 2006), which persists for the duration of the 
worms' lifespan. However, when Wolbachia is depleted from O. 
ochengi by oxytetracycline treatment, neutrophils are replaced by 
the classic helminthotoxic effector cell, eosinophils, which infiltrate 
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into the nodule and degranulate on the filarial cuticle, eventually 
resulting in worm death (Hansen et al. 2011). As this sequence of 
events is specifically triggered by disruption of the symbiosis rather 
than anthelminthic chemotherapy per se, wOo could be considered 
a defensive mutualist. This phenotype is now well recognized for 
Wolbachia in arthropods, and a precedent involving manipulation 
of the immune system of a third party exists for an insect pest of 
maize (Barr et al. 2010). Moreover, in L. sigmodontis, successful mi- 
gration of the infective larvae is also associated with immuno- 
modulation of mammalian effector cells by Wolbachia (Specht et al. 
2011). The abundance of TLR ligands with known or suspected 
chemo-attractant activity for neutrophils in the wOo proteome 
(WSP and PAL, and perhaps GroEL and SCOl) provides strong 
support for the hypothesis that Wolbachia not only stimulates the 
mammalian immune response but also can manipulate it in such 
a manner that the filarial host can survive for more than a decade. 

In conclusion, our global analysis of gene expression in wOo 
provides an unprecedented insight into the complexity of an or- 
ganism that can mediate different host interactions in separate 
tissue compartments, despite extensive genome reduction. Thus, 
we urge caution, following the recent explosion of knowledge re- 
garding Wolbachia phenotypes in arthropods, in any attempt to 
circumscribe the filarial- Wo/kzdz/a relationship in terms of a single 
process. In addition, we emphasize that our study was restricted 
to wOo within adult worms, and the stage-specificity of symbi- 
otic interactions has been recognized in other systems, such as 
the enterobacterial symbionts of entomopathogenic nematodes 
(Koppenhofer and Gaugler 2009). 

Methods 

Parasite material 

Nodules containing adult O. ochengi parasites were obtained from the 
skins of cattle slaughtered at Ngaoundere abattoir in the Adamawa 
Region of Cameroon (Wahl et al. 1994). Adult male and female 
worms were separated from excised nodules; the gonads were re- 
moved from the latter (selecting only nongravid individuals) and 
stored in PBS or RNA/ater (Sigma) at -80°C, whereas the males were 
preserved whole. The material was transported to the United King- 
dom on dry ice and stored at -80°C. 

Nucleic acid extractions 

O. ochengi genomic DNA was extracted using DNAzol (Invitrogen) 
from approximately 150 male worms in PBS or —20 mg female 
gonads in RNA/ater. For transcriptomic analyses, total RNA was 
purified from the worm material stored in RN Alater using TRI Re- 
agent (Sigma). For further details, see the Supplemental Methods. 

Genome sequence 

The wOo genome assembly was generated from the O. ochengi DNA 
extracts using an initial data set obtained on a Genome Sequencer 
FLX System (454 Life Sciences, Roche Diagnostics) and assembled 
in Newbler (454 Life Sciences), which was corrected for homo- 
polymer length errors with additional reads from a Genome Ana- 
lyzer IIx (Illumina) platform (final coverage, 70 X). The sequence 
was analyzed using Prodigal to call protein-coding genes, tRNAscan- 
SE (Lowe and Eddy 1997) to locate tRNA genes, MUMmer (Kurtz 
et al. 2004) to detect repeats, ISSaga (Varani et al. 2011) to iden- 
tify IS elements, BLAST (Altschul et al. 1990) and InterProScan 
(Zdobnov and Apweiler 2001) to assign putative functions to 
protein-coding genes, and ORTHOMCL (Li et al. 2003) to define 



orthologs shared with other Wolbachia strains. Metabolic path- 
ways were examined using the SEED (Overbeek et al. 2005) and 
Kyoto Encyclopaedia of Genes and Genomes (Kanehisa and Goto 
2000) databases. Pseudogenes were classified as regions with ho- 
mology with previously annotated bacterial genes but that showed 
evidence of disruption by premature stop codons or frameshift 
mutations. In addition to the automated correction of the 454 
assembly using Illumina genomic reads, putative pseudogenes 
were checked for data inconsistencies manually using both the 
Illumina data set and SOLiD transcript reads (see below). Finally, 
mass spectrometric data from protein samples (see below) were 
searched against six-frame translations of highly transcribed 
pseudogene fragments and expressed intergenic regions, as iden- 
tified by Cufflinks (Trapnell et al. 2010). For further details, see the 
Supplemental Methods. 

Synthesis and purification of cDNA 

Contaminating genomic DNA in the RNA samples was digested 
using a TURBO DNA-free Kit (Ambion) according to the manufac- 
turer's instructions. For first-strand cDNA synthesis, —2.5 juug RNA 
was reverse-transcribed using a Moloney murine leukemia virus 
reverse-transcriptase (RT) kit (Bioline), as specified by the manu- 
facturer. Second-strand synthesis was performed with E. coli DNA 
ligase, DNA polymerase I, RNase H, and T4 DNA polymerase in 1 X 
Second-Strand Reaction Buffer (all supplied by Invitrogen). The 
double-stranded (ds) cDNA was purified by phenol-chloroform 
extraction, dissolved in nuclease-free water, and stored at -20°C. 
For further details, see the Supplemental Methods. 

Transcriptome sequencing and analysis 

The ds cDNA samples (two biological replicates from male or fe- 
male worms) were used to generate standard SOLiD fragment li- 
braries and run as a multiplex on three SOLiD (version 4; Applied 
Biosystems) slide quadrants in total. The resulting reads where 
mapped to the finished genome using Bowtie (Langmead et al. 
2009). Transcript expression was analyzed using Cufflinks (Trapnell 
et al. 2010) to find novel expressed regions in the genome, whereas 
differential expression was calculated using reads counts generated 
from bam2rpkm freeware (obtained from SourceForge) and analyzed 
using edgeR (Bioconductor), with a negative binomial distribution 
model. To identify sRNA genes, noncoding regions exhibiting high 
levels of transcription were analyzed using Rfam (Griffiths-Jones et al. 
2003). For further details, see the Supplemental Methods. 

Sample preparation for proteomic analysis 

To provide Wolbachia-emiched samples for proteomic analysis, 
600 mg of fresh adult female O. ochengi was processed using the 
FOCUS Mitochondria kit (G-Biosciences) according to the manu- 
facturer's instructions for soft tissues. Proteins were separated un- 
der reducing conditions on a 15% polyacrylamide gel and visual- 
ized using a Colloidal Blue Staining Kit (Invitrogen). Thirty gel 
slices from a single lane (30 |xg total protein) were digested with 
trypsin, and five pooled fractions were analyzed by mass spec- 
trometry (MS). For further details, see the Supplemental Methods. 

Liquid chromatography-MS and data analysis 

Tryptic peptides were analyzed by high-resolution liquid chro- 
matography (LC) MS/MS on an LTQ-Orbitrap Velos (Thermo 
Fisher Scientific). Peptide identifications were established using 
the Mascot (Matrix Science) search engine against a custom data- 
base. Search parameters included a precursor mass tolerance of 10 
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ppm and fragment mass tolerance of 0.5 Da, while one missed 
tryptic cleavage was permitted. Carbamidomethylation was set as 
a fixed modification, and oxidation (M) was included as a variable 
modification. A significance threshold of P < 0.01 was applied to 
the Mascot ion score, and the FDRs were set at < 1%. For further 
details, see the Supplemental Methods. 

Data access 

The genome data have been submitted to the NCBI GenBank (http:// 
www.ncbi.nlm.nih.gov/genbank/) (accession PRJEA81837) and 
EMBL-Bank (http://www.ebi.ac.uk/embl/) (accession HE660029). 
The RNA-seq data are available from the Sequence Read Archive (ac- 
cession SRP014524), and the mass spectrometric data can be down- 
loaded from the Tranche repository at ProteomeCommons.org (hash 
+7EETqKbefFAa9amwlvZmw2VD36SYyrAjtRflh7fXkaXQeWTVpO/ 
HX3ihjTt21hoQrpqjVXxRon5f5TK74c6BU3TDXMAAAAAAAAElw==). 
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